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A coordination-based approach to elasticity of floppy and stiff random networks 
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We study the role of connectivity on the linear and nonlinear elastic behavior of amorphous 
systems using a two-dimensional random network of harmonic springs as a model system. A natural 
characterization of these systems arises in terms of the network coordination relative to that of 
an isostatic network 8z; a floppy network has 5z < 0, while a stiff network has Sz > 0. Under the 
influence of an externally applied load we observe that the response of both floppy and rigid network 
are controlled by the same critical point, corresponding to the onset of rigidity. We use numerical 
simulations to compute the exponents which characterize the shear modulus, the amplitude of non- 
affine displacements, and the network stiffening as a function of 5z, derive these theoretically and 
make predictions for the mechanical response of glasses and fibrous networks. 

PACS numbers: 62.20.de 



The mechanics of crystalline solids is a fairly well un- 
derstood subject owing to the simplicity of the underlying 
lattice which is periodic. In contrast, an understanding of 
the mechanics of amorphous solids is complicated by the 
presence of quenched disorder, often on multiple scales. 
Two structural properties affecting the elasticity of dis- 
ordered solids are their coordination, and the presence of 
different types of interactions between the constituents 
of vastly dissimilar strengths. In the case of weakly- 
coordinated covalent glass such as amorphous selenium, 
the backbone is floppy, i.e. it is continuously deformable 
with almost no energy cost, but weak interactions such 
as van der Waals are responsible for the non-vanishing 
elastic moduli. On the other hand, highly-coordinated 
covalent glasses such as silica, or amorphous particle as- 
semblies where the main interaction is radial, such as 
emulsions, metallic glasses or granular matter, the back- 
bone is stiff. In foams and fibrous networks which are 
made of low-dimensional structures such as filaments and 
membranes, there is a wide separation of energetic scales 
between stretching and bending modes. This leads to a 
range of curious mechanical responses in these materials 
including strongly non-affme deformations [H, 0, H, El & 0] 
and elastic moduli that can be quite sensitive to the ap- 
plied stress (uja. Despite several theoretical advances 
H, H, Qil HHnU ) a unified descriptions of these behav- 
iors remains to be given. Here we study the mechanical 
response of simple floppy and rigid systems as the coordi- 
nation is continuously varied and propose such a unifying 
approach. 

We start by recalling Maxwell's criterion for rigidity 
in a central force network [l3| by considering a set of N 
points in d dimensions, subject to iV c constraints in the 
form of bonds that connect these points. This network 
has Nd — N c effective degrees of freedom (ignoring the 
d(d + l)/2 rigid motions of the entire system), and an 
average coordination number z = 2N C /N. The system 
is said to be isostatic when the system is just rigid, i.e. 
the number of constraints and the number of degrees of 
freedom are just balanced, so that Nd = N c , and z = 2d. 
When z < z c , the network exhibits collective degrees of 



FIG. 1: (a) Stretching a zig-zag chain of springs that are freely 
hinged costs no energy until they are aligned; (b) A random 
network of such springs, when subjected to shear, exhibits 
a nonaffine displacement field (see text for definition). Here 
we show the result for 10000 random springs with coordina- 
tion z — z c — 4.0, subject to a shear strain 7 = 0.005; (c) 
Zooming in shows the presence of large correlated rotational 
deformations. 



freedom with no restoring force; these solutions are called 
soft modes. For such a network made of harmonic springs 
of stiffness k, the energy can be written as 
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where ftij is the unit vector going from i to j, and 5R{ 
is the displacement of particles i. Soft modes satisfy 
SE — 0, or equivalently (SRi ~ SRj) ■ — Vij. In 
Fig|TJa), we show the example of a one dimensional zig- 
zag structure that straightens without an energy cost un- 
til the external load does not couple to the soft modes, 
i.e. when the latter correspond to node displacements 
that are orthogonal to the load direction; the system then 
becomes stiff. 

To understand how this simple idea carries over to a 
floppy network, we created disordered two-dimensional 
networks with up to 10000 particles using biperiodic 
jammed configurations of bi-disperse particles [14]. Lin- 
ear springs were then used to connect neighboring par- 
ticles; starting with z sa 5.5 > z c , a family of net- 
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FIG. 2: (a) Dimensionless stress-strain curves for z = 3.8, 4.0, 
and 4.2. (b) Critical strain 7* vs. Sz = z — 4. Each point av- 
erages over 2 configurations, (c) Log-log plot of dimensionless 
stress-strain curve for 2 = 4. (d) Rescaled stress-strain curves 
for z £ [3.0, 4.99] shows that both floppy and stiff networks 
can be described in terms of the relative coordination Sz. 

works with lower coordination is then generated by pro- 
gressively removing links from the most connected pairs 
of nodes leading to isotropic networks with low hetero- 
geneity in both density and coordination, with a range 
of spring rest lengths with the ratios of extremes vary- 
ing by a factor of less than 2. The size of the small- 
est springs is I, and their stiffness k. These networks 
are different from those studied in rigidity percolation 
Q or self-organized networks @, where the fluctuations 
in coordination are dominant and can rigidify the sys- 
tem even if z < z c . In addition, to model the pres- 
ence of weak interactions, we introduce weak springs 
(of dimensionless stiffness k w = k w /k = 1CP 5 excepted 
when stated otherwise. ~ designates dimensional quanti- 
ties) with a number density p w which stabilize the sys- 
tem. We impose a pure shear deformation on the net- 
work incrementally, and minimize the system energy via 
a damped molecular dynamics method. The spatially av- 
eraged dimensionless Cauchy stress is then calculated as 
a" = a/k = Tlij fij ~ Ri)> where the sum is on 

all springs ij, V = V/l 2 is the dimensionless area of the 
simulation cell, = fij/kl is the dimensionless force 

vector in the link ij, and Ri = R/l is the dimensionless 
position of particle i. Fig|TJb,c) shows the response to 
shear for z = 4 and 7 = 0.05. 

In Fig.([2la) we quantify the elastic properties as a func- 
tion of parameters, and show the existence of three qual- 
itatively different stress-strain relations. For floppy net- 
works with z — z c = Sz < 0, there is a critical strain 7* 
separating a zero stress plateau and a strain-stiffening 
regime; this critical strain 7* is a function of the deficit 
in coordination number Sz and follows the scaling law 



7* ~ I Sz \^ , with (3 — 1, over nearly two decades up to 
z = 3 and a strain as large as 40%. For an isostatic sys- 
tem with Sz = 0, the system resists shear deformations 
nonlinearly as soon as 7 ^ 0; indeed the stress-strain 
relation is parabolic a ~ j 2 , as shown in Fig.([2]c). Fi- 
nally, for a rigid system with Sz > 0, a linear stress-strain 
regime can be identified; however both the value of the 
dimensionless shear modulus G = G/k and the extent 
of the linear regime vanishes as Sz — > 0, as we discuss 
quantitatively below. 

Since the scaling relations in Fig. (J5Jb,c) are reminis- 
cent of a critical point, we look for scaling functions on 
which all the stress-strain curves must collapse after a 
suitable re-scaling of the axis. Postulating G ~ Sz for 
z > z c , we write: 

a = G(Sz) 1 f ± ( 1 / 1 *)^Sz e 1 f ± ( J L ) (2 ) 

where the functions /+ (/_ ) characterize the stress-strain 
relation of rigid (floppy) networks when Sz > (Sz < 0). 
When the argument of / vanishes, the existence of both 
a linear and floppy regime implies that /+ (2) — > c, where 
c is a constant of order one, and f—(x) — ► 0. For strains 
much larger than the characteristic strain 7* , the shear 
modulus must depend on the strain but not depend sig- 
nificantly on the value of Sz. This implies the existence 
of a power-law form f(x) ~ x x , with y = 6/f3. Previous 
numerical and empirical studies [1,0, [lj| show that 9 = 1 
in two and three dimensions, and were justified theoreti- 
cally [13] ■ Our numerical results in Fig.([2]b) imply that 
f3 = 1 and Fig. (0c) implies \ = 1 , in agreement with the 
relation derived. To directly check the validity of Eq. 
we rescale the axes of Fig.([2la), and show the results in 
Fig.([2]d) vindicating our choice of scaling functions. 

Moving beyond the scaling properties of the stress- 
strain curves, we now consider the structure of the dis- 
placement fields shown in Fig.{T]b). The non-affine dis- 
placement field {SRn a (^)}i=i...N is defined as SR l n a = 
R l — R-a where R l is the final equilibrated configuration, 
and R l a is the affine configuration obtained after the ho- 
mogeneous strain 7 is imposed. We focus on the ampli- 
tude of the relative displacement among nearest neigh- 
bors, defined as SR^ a = (\SR\ a — SR J n a \) where the 
average considers all pairs of particles in contact. This is 
a measure of the rotation of individual springs in the non- 
affine field, and the variation in the relative displacement 
per unit strain defined as SV^~ a = dSR^ a (^)/dj charac- 
terizes this rotation as a function of strain increment. In 
FigJSa), we show that 5V£ a . ~ \Sz\^ with <j) = -1/2, i.e. 
as the rigid sub-skeleton or backbone becomes strongly 
constrained, larger and larger displacements are required 
to accommodate a given shear. The rigid backbone acts 
as a lever, whose amplification factor diverges as one ap- 
proaches the onset of rigidity corresponding to Sz = 0. 
This scaling has already been observed in assemblies of 
elastic particles above the rigidity onset 0]. 

A small affine deformation of amplitude 7 causes di- 
mensionless forces \SF) = {F 1 } on each node. Assum- 
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FIG. 3: (a) The relative non-affine displacement per unit 
strain as denned in the text 5V n . a . for 7 = 7*/10 and (b) 
the dimensionless shear modulus G vs. — Sz, in the floppy 
regime, for k w = 10 -J . The variation in (c) the dimensionless 
shear modulus G and (d) the velocity 8V^ a vs. Sz € [—1, 1], 
with kw = 1/300, p w = 4.1. We observe in particular the large 
particle velocity in the vicinity of the isostatic point 8z — 0. 



ing for convenience of notation that the stiffnesses and 
lengths of all spring are equal leads to F l = ' 
r • riij)riij, where T is the corresponding strain ten- 
sor, and the sum is on all the neighbor j of i. For a 
rigid system with Sz > 0, the non-affine displacements 
\SR) = {SR l n a }, which corresponds to the displacement 
of the particle along the direction of unbalanced forces, 
is: 



\8R) = M- l \5F) = ^\{8R{u)\5F)\5R{uj)) 



(3) 



where M is the dynamical matrix and \5R(u>)) is the nor- 
malized normal mode of stiffness to 2 . Then the relative 
displacement between neighboring particles is: 



SRij = SRi - 8R 3 = ^(5R{uj)\5F)SR t] ( 



(4) 



Following the justification of [11[, we shall assume that 
the contribution of the different modes are independent: 



J2±{5R(u)\SF) 2 {\\6R k 



(5) 



For weakly coordinated systems, we shall use the results 
of [HI which show that above some frequency uj* ~ Sz, 
(i) the density of states D{w) does not depend signifi- 
cantly on uj and (ii) above uj* , the modes are very het- 
erogeneous, so that the correlations of the displacements 
of the particles are weak, see also [l6| . For such modes 



(\\5R i3 (uj)\\ 2 ) ~ (||^-(a;)|| 2 ) = 1/N, the latter equality 
stemming from the normalization of the modes. To esti- 
mate (5R{uj)\SF) 2 = V^r)//, ')//, ) ■ - \ ■ ii,, ) 1 . 

we use the weak spatial correlation of the modes and 
treat the different terms as independent. Using ((Hij ■ F ■ 
ftij)] ) ~ 7 2 then leads to 

(SR(co)\SF) 2 ~ 7 2 J^[(SRi(co) - SR^lj)) • n %] )] 2 ~ 7 V 

ij 

(6) 

where we used the definition uj 2 — 1/2 [(8Ri(uj) — 
SRj(ui)) ■ riij)] 2 . Finally, in the large N limit we have 
XL 1 / N J dujD(uj), so that Eq.® yields: 
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UJ* Sz 

(7) 

leading to the relation SR^ a ~ H^-R^H ~ j/y/Sz. 

From this we may deduce the strain at which stiff- 
ening occurs. The small applied affinc deformation 7 
causes forces F % ~ 7, which lead to small non-affine dis- 
placements SR l n a » r yu l , where u 1 = lim 7 ^o^-R^. a ./7- 
Since the linear approximation SR l n a = ju 1 is not exact, 
there are small residual forces on the nodes, but an it- 
erative perturbative procedure can be used to determine 
the correction to the leading order result. These residual 
forces may be estimated following Pythagoras' theorem: 
the transverse relative displacement at a contact cause 
a strain and therefore a residual force of the order of 
^n. a 2 {l) ~ ^i 2 Sz 2 ^ . When this quantity becomes larger 
that F l *~ 7, the linear approximation breaks down. This 
occurs for some 7* ~ 5z@ ~ Sz~ 2 ^, yielding the relation 
(3 = —2cf)= 1 observed in our numerical simulations. The 
resulting divergence in the particle velocity explains why 
non-linearities occur at a very small strain close to the 
rigidity onset. 

The previous argument also yields the scaling form for 
the dimensionless shear stress G in the floppy regime. In- 
deed when a shear strain 7 is imposed, each weak spring 

stores a dimensionless energy of order 5E ~ k w SV^~ a , 
leading to G ~ p w k w SE/j 2 ~ p w k w /\Sz\ in the floppy 
regime, as observed in Fig([3]b). Obviously, this scal- 
ing is expected to fail near the rigidity threshold, and 
we expect a cross-over to occur when G ~ Sz as ex- 
pected for the backbone from our earlier scaling argu- 
ments. This allows us to define a characteristic coordi- 
nation scale u* r~j \J p w k w so that there is an associated 
critical strain 7* ~ u* ~ sj p w k w where our two esti- 
mates are of the same order. This defines three regimes 
for the mechanical response as a function of the relative 
coordination: (i) For Sz << —u*, G ~ (p w k w )/\Sz\ and 
7* ~ — Sz. The rigid backbone is not significantly de- 
formed, (ii) For —u* << Sz « u* , the energy is shared 
between strong and weak interactions, G ~ u* ~ %/ p w k w 
and 7* ~ u* . (iii) For Sz » u* , most of the energy is 
condensed in the rigid backbone, G ~ Sz and 7* ~ Sz. 
Fig([3Jc) shows the behavior of the shear modulus as the 
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coordination increases, and the cross-over from the floppy 
to the rigid network. To confirm this description, we 
compute SV^ a for different coordinations. The result is 
plotted in Fig([3Jd) for k w = 1/300. Although the ratio 
of stiffnesses between weak and strong springs is large, 
the intermediate region is of significant amplitude and 
vanishes only as the square root of this ratio, due to the 
"lever" effect induced by the backbone. 

We conclude with a brief discussion of our results and 
its implications. Using numerical simulations of a weakly 
coordinated network as an exploratory tool, we have 
shown that (i) rigid and floppy networks are controlled by 
the same critical point (ii) two exponents, 4> and 9 which 
characterize the amplitude of non-affine displacements 
and the shear modulus respectively completely character- 
ize the system; all other exponents describing the effects 
of non-linearities and the stiffness induced by weak inter- 
actions follow from these. Furthermore, near the rigidity 
threshold, the amplitude of the non-affine displacement 
rapidly increases and the material response is then char- 
acterized by a point on a two-dimensional phase diagram 
(Sz, u*). In the context of glasses, our model may de- 
scribe tetrahedral network glasses, such as silica, amor- 
phous silicon or water, where fluctuations of coordination 
are rare at reasonable pressures. Such networks, where 
the joints linking tetrahedra are soft, are marginally rigid 
(T7I |22| . and the dominant weak interaction rigidifying 
the system is the energy required to change the angle be- 
tween two adjacent tetrahedra, whose strength is small 
for silica, and much larger for silicon. Thus silica is ex- 
pected to behave effectively as a weakly-coordinated net- 
work, and exhibit large non-affine displacements, while 
amorphous silicon should not. In contrast, for chalco- 
genide glasses where the composition of atoms of differ- 
ent valence can be modified to control rigidity [18j, the 
precise topology of the covalent network near the rigid- 
ity threshold is still unsettled [2(3] ■ Our model is 



presumably the simplest approximation of such networks 
and gives a plausible explanation for the observation of 
a smooth cross-over in the elastic moduli near z c [23l |: 
although Van der Waals interactions are two orders of 
magnitude softer than covalent bonds, their effect is not 
negligible near z c due to the lever effect of the backbone 
that we have described, and further predicts that the 
non-affine displacement will be maximum near the rigid- 
ity threshold. In solid foams or stiff fiber networks, the 
different interactions at play (bending, stretching, cross- 
link rigidity...) cause a change in the effective coordi- 
nation. Since bending is a softer mode than stretching, 
those fibrous systems will be in general floppy. Our work 
yields the simple prediction that the strain 7* at which 
the system begins to stiffen and the amplitude of the non 
affine displacement will be anti-correlated, and that 7* 
will decrease as cross-links are added or fiber length is 
increased. At a quantitative level, we expect the corre- 
sponding exponents to depend on t he p articular struc- 
tural properties of the network [T(| [Hf. In particular, 
(i) if long fibers are present, non-affine displacements are 
enhanced, SV^.a. ~ I/I'M [l3|> an d non-linearities are 
expected to occur at a smaller strain, and (ii) if the link 
size is widely distributed, small links will effectively act 
as points of higher coordination, which will also tend to 
reduce 7*. More generally, our study suggests that the 
amplitude of the non-affinc displacements allows us to 
classify amorphous solids according to their closeness to 
criticality in the coordination number. 
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